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Abstract 

> 

CN , The high temperature asymptotics of thermodynamic functions of electro- 

magnetic field subjected to boundary conditions with spherical and cylindrical 

f^ ■ symmetries are constructed by making use of a general expansion in terms 

of heat kernel coefficients and the related determinant. For this, some new 
heat kernel coefficients and determinants had to be calculated for the bound- 
ary conditions under consideration. The obtained results reproduce all the 
Q^l asymptotics derived by other methods in the problems at hand and involve a 

few new terms in the high temperature expansions. An obvious merit of this 
approach is its universality and applicability to any boundary value problem 
correctly formulated. 
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I. INTRODUCTION 

The Casimir effect is one of the most interesting phenomena in quantum field theory. 
Since its discovery more than 50 years ago it attracted much attention. In the past years 
the interest intensified after its experimental verification reached the one percent level of 
precision ||l]-0]. 

The influence of temperature on the Casimir effect was an important topic since its first 
experimental demonstration |^ which had been done at room temperature. It was first 
shown in Ref. that the temperature influence was just below what had been measured. 
It is expected that the temperature contributions will be seen in the upcoming series of 
experiments. 

In quantum field theory, finite temperature effects can be described at equilibrium in the 
Matsubara formalism by imposing periodic (resp. antiperiodic for fermions) boundary con- 
ditions in the imaginary time coordinate. Technically this is very similar to the calculation 
of the Casimir effect for plane boundaries and can mathematically be described by the same 
Riemann and Hurwitz zeta functions. Another formula important for applications is the well 
known Lifshitz formula describing the interaction between plane dielectric bodies at finite 
temperature. In the case of nonflat boundaries the situation is, however, more complicated. 
In order to obtain the Casimir effect at finite temperature one has to know it at least at 
zero temperature, i.e., one has to know the spectrum of the corresponding operator. Even 
then complicated calculations are usually necessary and explicit results are rare. 

An opposite situation occurs with the asymptotic expansion of the Casimir energy at 
high temperature. It turns out to be determined to a large extent by local quantities which 
are much easier to obtain. In the paper 0, which did not receive the due attention, it was 
shown that this expansion can be written down in terms of the heat kernel coefficients and 
the functional determinant of the operator corresponding to the spatial part of the problem. 
For example, this method has been applied in 0-^ for the effective potential in curved 
spacetimes and for the Casimir effect for hypercuboids. 

The heat kernel coefficients for differential operators on manifolds with and without 
boundary are known to depend on the properties of this manifold only locally. This means 
that they can be represented as integrals over the manifold and over its boundary |T0[ 



whereby the characteristics like curvature enter as local functions. The calculation of these 
coefficients is a topic of its own and much progress had been achieved especially during the 
past decade (see the book [|l^, for example). Less is known about the determinant. It is not 
ascertain in general whether it is a local quantity. However, in several examples it is shown 
to be calculable much easier than the corresponding Casimir energy at zero temperature. 
As a consequence, the high temperature expansion of the Casimir energy can be calculated 



quite easily. This was emphasized in the recent review [0 where the basic formula (2.10) 
is taken from. 

In the present paper we apply these general formulas to some specific examples. First 
we consider parallel plates in order to demonstrate the technical tools on a simple problem. 
Then we consider the conducting spherical and cylindrical shells obtaining new terms in 
the asymptotic expansion. Eventually we consider the dielectric ball and cylinder where we 
restrict ourselves to the dilute approximation. These examples demonstrate the effectiveness 
of the method. 



An interesting application of the general formula is the discussion of the so called classical 
limit which had recently been considered in Ref. [0. It is understood as to take place if the 
internal energy which is connected with the free energy by Eq. (2.11) (see below the next 
section) tends to zero for T — ^ oo. This happens if the heat kernel coefficients with number 



n < ^ vanish. 



The first calculation of the leading contributions to the high temperature asymptotics of 
the Casimir energy for curved boundaries was given in Ref. |]I?| using the multiple refiection 



expansion. As it was noticed in the recent paper |]T5| the multiple refiection expansion can 
be used for the calculation of the heat kernel coefficients demonstrating the equivalence of 
both approaches up to the question of the determinant. 

The situation is to some extent different for boundaries with edges and corners. Here 
the application of Riemann and Hurwitz zeta functions seems to be more appropriate. A 



first example of this kind was given in Ref. [T^. The appropriate more general methods can 



be expected to be those given in Ref. [1/ 



The layout of the paper is as follows. In Sec. II the derivation of the high temperature 
expansions in terms of the heat kernel coefficients is briefiy given. In Sec. Ill the original 
setting of the Casimir effect, i.e. parallel perfectly conducting plates in vacuum, is considered 
and the high temperature asymptotics of the thermodynamic functions are derived in terms 
of the relevant heat kernel coefficients. In Sec. IV the high temperature asymptotics for 
electromagnetic field with boundary conditions on a sphere are obtained. In Sec. V the 
high temperature expansions are constructed for the boundary conditions defined on the 
lateral of a circular infinite cylinder. The heat kernel coefficients needed are calculated by 
making use of the respective zeta functions that have been obtained in an explicit form in 



terms of the Riemann zeta function in Ref. |18] and also by applying the results of Ref. [19 



The functional determinants entering the asymptotic expansions at hand are calculated by 



making use of the technique developed in Ref. pO[. The results obtained are compared with 



the high temperature asymptotics which have been derived for boundary conditions under 
consideration by other methods. The possible extension of the approach is discussed in the 
Conclusions (Sec. VI). 

The mathematical details of the calculation of the zeta determinants are presented in 
Appendix A for electromagnetic field subjected to boundary conditions given on a sphere 
and in Appendix B for the boundary conditions defined on the lateral of an infinite circular 
cylinder. 

II. HEAT KERNEL COEFFICIENTS AND HIGH TEMPERATURE 

EXPANSIONS 



Let the dynamics of quantum field be defined by the operator 



1 d^ 

where A is not of necessity the Laplace operator, but an elliptic differential operator de- 
pending only on space coordinates. The free energy F of the field is determined by the zeta 
function Ct("S) corresponding to the Euclidean version of the operator ( p^.lj ) 



^ = -|Ct(0). (2.2) 

Here T is the temperature measured in energy units (the Boltzmann constant k^ is assumed 
to be equal to 1), and the zeta function Ct('S) is defined in a standard way 

Ct(^)= E EK + 4"'' (2.3) 

m=-oo {k} 

with Qrn = 2TTmT/h being the Matsubara frequencies and ujI/c^ standing for the eigenvalues 
of the operator —A in Eq. (|2.1| ) 



UJ 



Aifkix) = -iM^)- (2.4) 



c 



The characteristics of the quantum field system with dynamical operator (|2.1| ) at zero 



temperature are determined by the zeta function ({s) associated with the operator —A 

C(^) = E^."'^- (2-5) 

{k} 

From the mathematical point of view the zeta function ({s) corresponding to the space part 
of the operator ( p.lD is, undoubtedly, a simpler object than the complete zeta function Ct('5) 
because the definition (|2.3| ) involves an additional sum over the Matsubara frequencies. 
Here a natural question arises whether one can gain knowledge of the quantum field at 
nonzero temperature possessing only the zeta function ({s). In Ref. [^ it was shown that 
proceeding from the zeta function ({s) one can deduce the high temperature asymptotics of 
the thermodynamic functions such as Helmholtz free energy, internal energy, and entropy. 
Let us remind briefiy the derivation of these asymptotics. By making use of the formula 

\-' = —- / dtt'-^e-^' (2.6) 

r(s) io ^ ^ 

the zeta function (|2.3| ) can be represented in the form 



X roo 
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-*- {^) •'^ m=~OD {k} 



The term with m = in this formula gives the zeta function ( |2.5| ). In the remaining terms 
we substitute the heat kernel K(t) of the operator —A by its asymptotic expansion at small t 

m^T.e-''^j^, E a„r + .... (2.8) 

{A;} ^^'^'') n=0,l/2,... 

As a result we arrive at the following asymptotic representation for the complete zeta func- 
tion Ct(s) 

. . ^ .. ^ 2 ^ f h Y''^^^'"T(s- 3/2 + n) ^ , , , , 

CtW-CW+P^__J,^;4^J r(l) ' Cr(2. + ^"-3). (2-9) 



where Cr{s) is the Riemann zeta function. Taking the derivative of the right hand side of 
Eq. (|2.9|) at the point s = and substituting the result into Eq. ( p.2|) one obtains the high 
temperature expansion for the free energy 



F(T) 
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(2.10) 



Here 7 is the Euler constant. The argument of the logarithms in expansion ( p.lO|) are 
dimensional, but upon collecting similar terms with account for the logarithmic ones in C'(0) 
it is easy to see that finally the logarithm function has a dimensionless argument, at least for 
a2 = 0. Let us note that according to the definition ( |2.8| ) the heat kernel coefficients in our 
consideration are dimensional, because the frequencies Uk have the dimensionality [time]~^. 
The asymptotic expansions for the internal energy U{T) and the entropy S{T) are de- 
duced from Eq. ( |2.1CI| ) employing the thermodynamic relations 
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(2.11) 
(2.12) 



Substituting the expansion ( |2.10|) into Eqs. ( p.ll|) and ( p.l2| ) one arrives at the asymptotics 
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(2.13) 
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(2.14) 



In Eq. ( P^.13| ) the term proportional to 02 contains the logarithm of dimensional quantity: 
[h/T] = [tinie]"^. This is the result of the arbitrariness arising in from the ultraviolet 
divergences in the case of 02 7^ (see Ref. |21[] for a more detailed discussion). Unlike this 
situation, collecting the logarithm functions in the a3/2-term and in C{0) in Eq. ( p.l4| ) leads 
to a dimensionless argument of the logarithm in the final expression. 



It is worth noting that the zeta determinant of the operator —A (i. e. C'(0)) does not enter 
the asymptotic expansion for the internal energy ( 2.13| ). Therefore this high temperature 
expansion is completely defined only by the heat kernel coefficients. In view of this, the first 
term in the asymptotics of the free energy in Eq. ( p.lO|) is referred to as a pure entropic 
contribution. Its physical origin is till now not elucidated. 



III. PERFECTLY CONDUCTING PARALLEL PLATES IN VACUUM 



In this section we demonstrate the application of the high temperature expansions ( p. 101) , 
( p.l3| ), and ( |2.14| ) to a simple problem of electromagnetic field confined between two perfectly 
conducting parallel plates in vacuum. First, we briefiy recall how to construct the zeta 
function in this problem. 

As well known, for example, from the theory of waveguides and resonators [^ the vectors 
of electric and magnetic fields in the problem at hand are expressed in terms of the electric 
(n') and magnetic (11") Hertz vectors, each having only one nonzero component 11^ and 
n^ satisfying, respectively, Dirichlet and Neumann conditions on the internal surface of the 
plates. The functions 11^ and 11^ obey the equations 






+ v^ n 
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(3.1) 



where uj is the frequency of electromagnetic oscillations, V^ stands for the two-dimensional 
Laplace operator for the variables (x, y) = x. The separation of variables results in the 
following solution 



, ,'miz\ 
II^(x, z) = exp(«kx) sm ( I 
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n = l, 2, 



ri = 0, 1, 2, 



(3.2) 



where a is the distance between the plates. Hence, the states of electromagnetic field with 
the energy hun, n > 1, are doubly degenerate, while the state with the energy hujo = Tick 
is nondegenerate. 

With allowance for this the zeta function in the problem under consideration is given by 
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(3.3) 



where L^ and Ly are the dimensions of the plates. 

For a correct definition of the integral in this formula in the small k region the photon 
mass fi is introduced (infrared regularization) . At the final step of calculations one should 
put /i = 0. On integrating in Eq. (|3^ ) and substituting the sum over n by the Riemann 
zeta function one arrives at the result 



The zeta function 
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gives the well-known value for the Casimir energy 



(3.4) 



Er 



or for its density 
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where V^ = aL^ Ly. 



(3.5) 



(3.6) 



In order to construct the high temperature expansions (|2.1(]|) , (p.l3|) , and ( p.l4|) the heat 
kernel coefficients for the system under consideration should be obtained by making use of 



the zeta function (|3]J). 

The zeta function ( |3.4| ) or, in the general case, (p.5|) and the corresponding heat kernel 
"M are related via the Mellin transform 
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This enables one to express the heat kernel coefficients a^ in terms of the values of the zeta 
function at the respective points 
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(3.8) 



Substituting Eq. ( p.4|) into Eq. ( p.8|) we obtain for perfectly conducting parallel plates 
only one nonzero coefficient Qq 



flo 



4. 



(3.9) 



where V = L^Lya is the volume of the space bounded by the plates.Q This is just an 
illustration of the well-known fact that for flat manifolds without boundary or with flat 



boundary all the heat kernel coefficients except for oq vanish |^3[. It should be noted here 
that we are considering only electromagnetic field confined between the plates and do not 
take into account that field outside the plates. 

From Eqs. ( p.l3| ) and (|3.9|) it follows that the density of internal energy has the following 
high temperature asymptotics 
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(3.10) 



^For obtaining the vanishing a^^ ^ coefficient it is important to take into account the second term 
in Eq. (p.4|) which depends on the photon mass //. 



where a is the Stefan-Boltzmann constant 
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(3.11) 



Recall that in our formulae we put /cb = 1, that is, the temperature is measured in energy 
units. The transition to degrees is performed by the substitution T —^ k-^T . 

When calculating the high temperature asymptotics of the free energy (|2.10|) and the 
entropy ( p.l4| ) one needs to derive C'(0) for the zeta function ( p.4| ). Keeping in mind that 
Cr(— 2) = it is convenient to use here the Riemann reflection formula 



which yields 



From here we deduce 



2^-^ T{s) Ck{s) cos(7r s/2) = ti\^{1 - s) 



U2S-2) ^ -.^ + 0(.2). 



C'(0) = ^Cr(3) = -^Cr(3). 
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(3.12) 
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Insertion of Eqs. ( p.9| ) and p.l4| ) into Eq. ( |2.10| ) gives the following high temperature be- 
haviour for the density of free energy 
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(3.15) 



As was noted above, we are considering only electromagnetic field between the plates. There- 
fore when calculating the Casimir forces one should drop the last term in Eq. (|3.15|) since its 
contribution is canceled by the pressure of the black body radiation on the outward surfaces 
of the plates. As a result the high temperature asymptotics of the Casimir force, per unit 
surface area, attracting two perfectly conducting plates in vacuum is 



J" 



T 



Ana^ 



Cr(3). 



(3.16) 



Usually in the Casimir calculations the contribution of the free black body radiation is 



subtracted from the very beginning ||2J 



It is interesting to note that the Casimir force ( p. 16 ) and the first term on the right 
hand side of Eq. ( p.l5|) are pure classical quantities because they do not involve the Planck 
constant h. These classical asymptotics seem to be derivable without appealing to the notion 
of quantized electromagnetic field. The classical limit of the theory of the Casimir effect is 
discussed in a recent paper ]TB|] . 

Employing Eqs. ( 2.12|) and ( 3.15 ) one arrives at the high temperature behavior of the 
entropy density 



S{T) 
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Cr(3) 2T3^ 
87ra3 45c^h^ 



(3.17) 



It is worth noting that the corrections to Eqs. ( p.lO|) , ( |3.15| ), (|3.17|) are exponentially small. 



The example considered shows that the zeta function of the spatial part of evolution op- 
erator really enables one to obtain the high temperature asymptotics of the thermodynamic 
functions in a straightforward way. In the subsequent sections we shall consider quantum 
fields defined on manifolds with boundaries possessing spherical or cylindrical symmetries, 
when the relevant zeta functions cannot be obtained in a closed form. Furthermore in these 
cases the spectrum of the operator —A is not known explicitly. And nevertheless the method 
proposed is applicable to these cases also. 

IV. THERMODYNAMIC ASYMPTOTICS FOR ELECTROMAGNETIC FIELD 
WITH BOUNDARY CONDITIONS ON A SPHERE 

In the present section we consider electromagnetic field subjected to three types of bound- 
ary conditions on the surface of a sphere: i) an infinitely thin and perfectly conducting 
spherical shell; ii) the surface of a sphere delimits two material media with the same velocity 
of light; iii) a dielectric ball placed in unbounded dielectric medium. In order to obtain 
the heat kernel coefficients determining the high temperature asymptotics (p.lO|) , ( p. 131) , 



and (|2.14|) it is convenient to use the explicit representation of the relevant spectral zeta 



functions in terms of the Riemann zeta function. These formulae were derived in our recent 
paper [l^ by taking into account the first two terms of the uniform asymptotic expansion 
for the product of the modified Bessel functions Ivivz) K^^uz). 

A. Perfectly conducting spherical shell 

We take advantage of Eq. (2.26) in Ref. |T^ substituting there the variable s by 2s 
and recovering the explicit dependence on the velocity of light c. The latter results in the 
replacement of the sphere radius by R/c: 

1 fR\^' 



where 



({s) ^ 4 (^j ^ (1 + s) {2 + s){ (2^+2^ - l) Cr(1 + 2s) - 2i+2« 

+ q{s) [(2^+2^ - l) Cr(3 + 2s) - 23+2S] ^ _ _ _| ^ (4^^^ 



q(s) = (480 + 1736 + 2016 s^ + 568 s^), (4.2) 

yv ; 3g4Q V ;, v ; 



and R is the radius of a sphere. The terms omitted in Eq. (|]T) are of the form 

,2(k+s)+i _ 1 Wj^(2 A; + 2 s + 1) - 2^ (k+^)+^] , /t = 2, 3, 4, . . . , (4.3) 



qk{s) [(2 



where qk{s) stand for some polynomials in s. 

Analysis of Eqs. (|4.1| ) and (|4.2| ) shows that the zeta function (|4.1| ) for a perfectly conduct- 
ing spherical shell enables one to find the exact values of the first six heat kernel coefficients, 
namely: 

9 



7r3/2 c2 



ao = 0, ai/2 = 0, ai = 0, 03/2 = 2 7r^/^ aa = 0, 05/2 = — - -^. (4.4) 



Taking into account the structure of the omitted terms ( [4 .31) it is easy to see that 

a, = 0, j = 3,4,5,.... (4.5) 



Having obtained the heat kernel coefficients ( [4. 4)) and (|4.5|) we are in position to construct 
the high temperature asymptotics of the internal energy of electromagnetic field by making 
use of Eq. (^J^ ) 



U(T)^--(—]^^ + 0(T-^). (4.6) 

^ ' 4 \ R 1920 T ^ ' ^ ^ 



Applying the technique developed in Ref. |25| more terms to this expansion can be easily 
added. 



In order to write the asymptotic expansions ( p. 101 ) and ( |2.14| ) the derivative of the zeta 



function at the point s = should be calculated. Equation ( [4.1|) gives an approximate value 
for C'(0) 

('(0) = ^+ln2 + — CR(3)-- + - In - = 0.38265 + - In-. (4.7) 

^ ^ ^ 2 16 ^^ ^ 8 2 c 2 c ^ ^ 

The terms omitted in ( [4.1|) will render precise only the ffist term in the final form of this 
expression, while the second term (1/2) ln(i?/c) will not change. The exact value of C'(0) is 
calculated in Appendix A 

('(0) = -_I + iln2 + 6 Cr(-1) + -- + - In - + ln2 + ^ 
s V ; 2 2 6 ^^^ ^V82 c 2 

1 R 

= 0.38429 + -In-. (4.8) 

2 c 

It is worth noting that the expression in the round parentheses, being multiplied by ^^, is 
exactly the value of C'(0) for a compact ball with continuous velocity of light on its surface 
(see Eq. (|4.20|) in the next subsection). As a result we have the following high temperature 



asymptotics of the free energy and the entropy in the problem in question 

S(r). 0.44215 + i 1„ ff -5I5 (|^)%0(r-). (4.10) 

The expression ( [4.9[ ) exactly reproduces the asymptotics obtained in Ref. |T^ by making use 
of the multiple scattering technique (see Eq. (8.39) in that paper). We have not calculated 
the coefficient 07/2, therefore we do not know the sign of the T~'^-correction in ( [4.91) . In Ref. 



14 1 it is noted that this term is negative. 
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In Eqs. ( [4.6|) , ( [4.9| ), and ([4.10| ) the large expansion parameter is actually a dimensionless 



'temperature' r = RT/{hc). Therefore the same formulae describe the behavior of the 
thermodynamic functions when R ^ oo and temperature T is fixed. 

The high temperature asymptotics of the thermodynamic functions derived by making 
use of the general expansions ( |2.1CI| ), ( p.l3| ), and (|2.14| ) contain the terms independent of 



the Planck constant h or, in other words, classical contributions (see Eqs. ([4.6|) , (|4.9| ), and 
( [4.1(j| )). This is also true for the high temperature limit of the Casimir force calculated per 
unit area of a sphere 

^m~-^_5£m = ^__f^V^ ^ + 0(r-) (4 11) 

^ '" inIP dR 16irff \R) 4irflM920r ^ ' ^ ' 

The leading classical term in the asymptotics (|4.11|) describes the Casimir force that 
seeks to expand the sphere. The quantum correction in this formula stands for the Casimir 
pressure exerted on the sphere surface. 

In Eqs. (4^), (|4.9| ) and ( |4.10 ) the Stefan-Boltzmann terms proportional to T^ are absent 



because the contribution of the Minkowski space was subtracted from the very beginning in 



our calculations [|I^. As a result we obtain the vanishing heat kernel coefficient oq which, 
in general case, is equal to the volume of the system under study |T^. Therefore our 
results describe only the deviation from the Stefan-Boltzmann law caused by the perfectly 
conducting sphere. 

The vanishing of the coefficients ai/2 and ai in the problem at hand can be explained 
by taking into account the general properties of the heat kernel coefficients |T^ and by 



making use of the results obtained in Ref. [25|. As known p6| the solutions to the Maxwell 



equations with allowance for a perfectly conducting sphere are expressed in terms of the two 
scalar functions that satisfy the Laplace equation with the Dirichlet and Robin boundary 
conditions on internal and external surfaces of the sphere. In view of this one can write 



an = a^^ + a^_ + a^_^ + a^_, n = 1/2, 1, . . . , (4.12) 

where subscribe plus (minus) corresponds to internal (external) region and the rest notations 
are obvious. In Ref. [E3i it was found that 



a?/,^ = -2n^l^R^ = a?/,_, af^,^ = 2n^/'R' = a%_ 



a. 



D lSvtR ^r _ -p l67r-R lA-\'x\ 



As a result we have 



ai/2 = ai = 0. (4.14) 



Having calculated the corrections to the Stefan-Boltzmann law one should naturally 
discuss the possibility of their detection. The ratio of the leading term in Eq. ( |4.(j| ) to the 
internal energy of black body radiation in unbounded space given by the Stefan-Boltzmann 
law ( |3.10| ) is proportional to r~^. Already for r ~ 10 the corrections prove to be of order 
10~^. The same value of r can be reached by varying the scale of length R in the problem 
under consideration or by respective choice of the temperature T. Keeping in mind the 
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value of the conversion coefficient ch = 197.326 MeV fm = 0.229 K cm p^ w^ obtain the 
following estimations. For R ~ 10~^^ cm (a typical hadron size) the temperature T should 
satisfy the inequality T ^ 200 MeV in order to apply the asymptotics found. For i? ~ 1 cm 
we have T 3> 0.229 K and for i? ~ 7 • lO^'' cm (radius of the Sun) the range of applicability 
of the asymptotics at hand extends practically to any temperature value T ^ 10^^*^ K. Here 
we shall not go into the details of a concrete experimental equipment that enables one to 
observe the calculated corrections to the Stefan-Boltzmann law confining ourselves to the 
estimations presented above. 

B. Compact ball with equal velocities of light inside and outside 

Let us consider the spherical surface that delimits the media with 'relativistic invariant' 



characteristics, i.e., the velocity of light is the same inside and outside the sphere [Q. In 



this problem there naturally arises a dimensionless parameter [^ 

e = (^i^Y = (ti-i^)\ (4.15) 

where ei and 62 (yUi and /X2) are permittivities (permeabilities) inside and outside the sphere. 
As usual we perform the calculation in the ffist order of the expansion with respect to ^^. 

In order to derive the zeta function for the boundary conditions at hand one should 
multiply Eq. ( [4.1|) by ^^ and replace there q{s) by the polynomial 



pis) = -- 



l-^(3 + s) + ^(3 + .)(4 + .)-^(3 + .)(4 + .)(5 + .) 



(4.16) 

The zeta function, obtained in this way, affords the exact heat kernel coefficients up to 03 
ao = 0, ai/2 = 0, ai = 0, 03/2 = 2 tt^/^^^, 02 = 0, 

Q5/2 _ .2 c^ p(-i) _ n r4 1 71 



With allowance of the structure of the omitted terms in Eq. (|4.1| ) we can again deduce that 

a,=0, j = 3,4,5,.... (4.18) 

Substitution of these coefficients into Eq. ( p.l3| ) gives the following high temperature behav- 
ior of the internal energy in the problem under consideration 

U{T)c^^^j + 0{T'^). (4.19) 

The value of C'(0) is calculated in Appendix A 



C'(0) = e' (^- g + 2 In - + In 2 + ^ j = e' (^0.35676 + ^^^-)- (4-20) 

It is this value that is supplied by Eq. ( [4.1D after changes specified above and with allowance 
for that p(-l) = 0. 
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By making use of Eqs. ( |2.10|) , ( |4.17|) , and ( |4.2CI|) we deduce the high temperature asymp- 
totics for free energy 



FiT) = -ej{l + \nA-^-^ + ^^T\n^ + 0{T-% 
= -e |o.71352 + ^ T In ||; + 0{T-'). 
The entropy in the present case has the following high temperature behavior 

he 



(4.21) 



S{T) 



Ul+7 + ln4-- 



^ (1.71352 




+ 0{T- 



(4.22) 



The asymptotics ( [4.19|) and (|4.21| ) completely coincide with the analogous formulae obtained 
in Refs. [0,^] by the mode summation method combined with the addition theorem for 
the Bessel functions. 



In Ref. |3^ the exact expression has also been derived for the internal energy in the 
problem at hand (see Eq. (3.22) in that paper). This formula gives only exponentially 
suppressed corrections to the leading term (^4.19|) 



T r 

U{T) ~ i^- [l + 2(4t^ + 4t + l)e 



-At 



(4.23) 



where t = 2ttRT. We have used here the relation between ^^ and An: C,'^ = An'^/4 (see Eq. 
(3.12) in Ref. 0]). The asymptotics ( ^.23| ) implies in particular that in reality in Eq. ( |4.19| ) 
there are no corrections proportional to the inverse powers of the temperature T. From 
here it follows immediately that all the heat kernel coefficients with integer and half integer 
numbers equal or greater than 3 should vanish 



0, J = 3, 7/2, 4, 5/2, 5, 



(4.24) 



(compare with Eq. ( [4.18| )). In view of this the sign denoting the omitted terms in Eqs. 
(|1|), ( ^:2ll ), and (|t2^ ) should be substituted by Oie-^""^^). 



C. Dielectric ball in unbounded dielectric medium 



The zeta function for electromagnetic field in the background of a pure dielectric ball 
(/ii = /i2 = 1, £^1 7^ £2) has not been obtained in an explicit form. In Ref. [^] the heat 
kernel coefficients up to 02 in this problem were found. Here we use the results of this paper 
confining ourselves to the Arz^-approximation, where An = rii — n2 = nin2 (c2 — ci)/c ~ 
(c2— Ci)/c, Hi and q are the refractive index and the velocity of light inside {i = 1) and outside 
{i = 2) the ball, and c is the velocity of light in the vacuum: rii = y/el, Ci = c/rii, i = 1,2. 
It is assumed that Ci and C2 differ from c slightly, therefore C2—C1 and An are small quantities. 
In view of this we have 
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ao 

«l/2 = 

Oi ~ 0, 



WR 



.3C2 






Stt^ (An 



2An2 






c 



«3/2 



TT 



/ 2 2\2 

3/2i£lZ_£2j -- 7r3/2Ar)2 



02 ^ 0. 



(4.25) 



The coefficients ai and 02 equal zero only in the An^-approximation considered here. In the 
general case they contain terms proportional to An'^, where /c > 3. 

Allowance for one more term in the uniform asymptotic expansion of the modified Bessel 



functions, as compared with the calculations in Ref. |21[], gives the next heat kernel coefficient 

a^/2 25 c^ 



(47r)3/2 2688 R^ 



An^ 



(4.26) 



Correcting the mistake made in [^ we state that this coefficient has no contributions 
proportional to An^, and in the An^-approximation one has to put 



^5/2 



0. 



(4.27) 



Making use of the technique developed in Ref. |2^ one obtains the following expression 
for the derivative of the zeta function for a pure dielectric ball at the point s = (see 
Appendix A) 



C'(o) 



An^ 



7 , ^ , . 

--+ln — + ln4 + 7 

8 c 



(4.28) 



Before turning to the construction of the high temperature asymptotics in the problem at 
hand by making use of the general formulae (|2.10|) , ( |2.13|) , and (|2.14|) the following remark 
should be done. When considering the electromagnetic field in the background of a dielectric 
body in the formalism of quantum electrodynamics of continuous media, as a matter of 
fact one is dealing with a system consisting of two objects: electromagnetic field plus a 
continuous dielectric body. It is important that this body is described (phenomenologically) 
only by respective permittivity without introducing into the Hamiltonian special additional 
dynamical variables. As a result the zeta function and the relevant heat kernel coefficients 
calculated in this formalism also describe both electromagnetic field and dielectric body. 
When we are interested in the Casimir thermodynamic functions in such problems we have 
obviously to separate in the general expressions the contributions due to the dielectric body 
itself 0. 

Let us turn to such separation procedure in the high temperature asymptotics for a 



dielectric ball. Following the reasoning of Refs. ||3^ , |35| we divide the Helmholtz free energy 
of a material body with volume V and the surface area S into the parts 



F = V f + Sa + F^ 



Cas) 



(4.29) 



where / is the free energy of a unit volume of a ball, a denotes the surface tension, and 
Fcas is refered to as the Casimir free energy of electromagnetic field connected with this 
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body and having the temperature T. In this way we obtain the following high temperature 
behavior of the free energy F{T) in the problem at hand 



F(r)^ao . 



11^ 



J3 



where oq and ai/2 are defined in Eq. ( ^.25| ) and 

i^Cas(T) - 



CR(3) + Fca.(T), 



(4.30) 



An^ / 4Ti? 7\ ,^,^ 2x 



(4.31) 



The high temperature asymptotics for the Casimir internal energy and for the Casimir 
entropy can be derived by making use of the respective thermodynamical relations ( p^.ll ), 



Ar?^ 



-2\ 



^ ,^, An2 /I ^ ART 



0{T 



-3^ 



(4.32) 
(4.33) 



It is worth comparing these results with analogous asymptotics obtained by different 
methods. In Ref. |3^ at the beginning of calculations the first term of expansion of internal 
energy (|4.32| ) was derived. The subsequent integration of the thermodynamic relation (p.ll|) 
gave the correct coefficient of the logarithmic term in the asymptotics of free energy ( [4.311 ). 
In a very recent paper |^ Barton managed to deduce the asymptotics ( [4.31|) - ( [4.331) . 
One should keep in mind that our parameter An corresponds to 2'Kan in the notations of 



Ref. [35 



The asymptotics ( |4.31[ )-( ^.33[ ) contain the /^-independent terms. As far as we know the 
physical meaning of such terms remains unclear. 

Preliminary analysis of a complete expression for the internal energy of a dielectric ball 



(see Eqs. (3.20) and (3.31) in Ref. |30]) shows that probably there are only exponentially 
suppressed corrections to the leading term ( [4.32[ ). In that case in addition to Eq. ([4.27[) all the 
heat kernel coefficients with number greater than 3 should vanish in the An^-approximation. 



V. THERMODYNAMIC ASYMPTOTICS FOR ELECTROMAGNETIC FIELD 
WITH BOUNDARY CONDITIONS ON AN INFINITE CYLINDER 



The calculation of the vacuum energy of electromagnetic field with boundary conditions 
defined on a cylinder, to say nothing of the temperature corrections, turned out to be 
technically a more involved problem than the analogous one for a sphere. Therefore the 
Casimir problem for a cylinder has been considered only in a few papers [|T^ , p!8| , p!^ ,^6H39| . 
We again examine three cases: i) perfectly conducting cylindrical shell; ii) solid cylinder with 
ci = C2; iii) dielectric cylinder when ci 7^ C2. Here we shall use the results of our previous 
papers [p!^ , pJ|] . 
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A. Perfectly conducting cylindrical shell 

In Ref. [|^ the first two terms in tlie uniform asymptotic expansion of the product of the 
modified Bessel functions In{nx) Kn{nx) were taken into account. As a result the spectral 
zeta function in the problem under consideration was represented as an expansion in terms 
of the Riemann zeta functions CR(2(fc + s) + 1), A; = 0, 1, 2, ... . With allowance for the 
first two terms in this expansion the zeta function is given by 

as) = Z^{s) + Z2{s) + Z^{s). (5.1) 

Here the function Zi{s) stands for the contribution of zero orbital momentum with proper 
subtraction 

1 

/^n(y) = y {In{y) Kn{y))' , t{y) 



The function Z2{s) is generated by the first term of the uniform asymptotic expansion 



R2«-i rn /2 -I- si 

.-Tr-^s (1 - 2^) (3 - 2^) [2 Cr(2. + 1) + 1] ^^7 "^ 
GAJttc'^'^ T{s] 



Z2{s) = ,, ;^,2s (1 - 2^) (3 - 2s) [2 Cr(2s + 1) + 1] ^;,„, ' . (5.3) 



The function Z^ corresponds to the second term of the uniform asymptotic expansion 

^3(^) = ,.,,n r ^1 - 2s) (3 - 2.) (784 s' - 104 s - 235) ^^ ^ ^ U2s + 3) . (5.4) 
61440 VTT T[s) 

The function Zi{s) is defined in the strip —3/2 < 3fJs < 1/2, while the functions Z2{s) 
and Z^{s) are analytic functions in the whole complex plane s except for the points, where 
r(s) and C,k[s) have simple poles. In order to find the heat kernel coefficients Oq, ai/2, and 
Oi trough the relation ( p.8|) one needs the zeta function defined in the region 1/2 + e < 
3? s < 3/2 + e with e being a positive infinitesimal. However in this region Eq. ( p.2[ ) is not 
applicable directly due to the bad behaviour of the integral at the upper limit. In the most 
simple way we can overcome this difficulty as in the case of perfectly conducting plates by 
introducing the photon mass /i at the very beginning of the calculation and making then 
the analytic continuation of the zeta function to the points s = 1/2, 1, 3/2. Upon taking 
the residua at these points one should put /i = 0. 

With regard to all this and using the relation (p.8|) we find the heat kernel coefficients 



ao = 0, ai/2 = 0, fli = 0, 02 = 0. (5.5) 

The vanishing heat kernel coefficient 02 implies that the zeta regularization gives a finite 



value for the vacuum energy in the problem at hand |[18|j37[] . The coefficient 03/2 is determined 
by the function ^'2(5) only (see Eq. (p.3|)) 



"^/^ - ' (5.6) 



(47r)3/2 64 i? 
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The coefficient 05/2 is defined by the function ^3(5) given in Eq. ( |5.4| 

0-5/2 _ 153 c^ 
(47r)3/2 ~ 8192^' 



(5.7) 



The calculation of the next heat kernel coefficients 03, 07/2, . . . would demand a knowledge 
of the additional terms in the expansion of the spectral zeta function in the problem under 
consideration in terms of the Riemann zeta function. These terms are proportional to 
Cr(2/c + 2s + 1) with /c = 2, 3, . . . , and may be obtained employing the technique developed 
in Ref . |TB[ . Analyzing the position of poles for these Riemann zeta functions it is easy to 
show that, as well as in the spherical case, we have 

aj = 0, j = 3, 4, 5 . . . . 

The zeta determinant entering the high temperature asymptotics of free energy (|2.10|) 
and entropy ( p. 14 ) is calculated in Appendix B 



Now we are able to construct the high temperature expansions of the thermodynamic 
functions in the problem under consideration. For the free energy we have 

T 3T , RT 51 Ti^c^ ,^, ,, , , 

FiT) . -0.22924 --_,„_- ^^ _ + 0(T-'). (5.9) 



When comparing Eq. (|5.9| ) with results of other authors one should remember that all the 
thermodynamic quantities that we obtained in this section are related to a cylinder of unit 
length. The high temperature asymptotics of the electromagnetic free energy in presence of 



perfectly conducting cylindrical shell was investigated in Ref. [|14[. To make the comparison 
convenient we rewrite their result as follows 

F(r). -0.10362 1 -|^lu|I. (5.10) 



The discrepancy between the terms linear in T in Eqs. ( ^.9|) and (|5.10|) is due to the double 



scattering approximation used in Ref. [T^ (see also the next subsection). Our approach 
provides an opportunity to calculate the exact value of this term (see Eq. (|5.9|) ). 
And finally, making use of the general formulae ( |2.13| ) and ( p.l4| ) we derive 



S(T).^lI51! + ^l„^_^5^4^ + 0(r-). ,5.12) 

^ ' R 64i? 2hc 196608 R^T^ ^ ' ^ ^ 
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B. Compact cylinder with ci = C2 and with ci 7^ C2 



Here we consider the boundary conditions for electromagnetic field of two types: i) a 
compact infinite cylinder with uniform velocity of light on its lateral surface, ii) a pure 
dielectric cylinder with ci 7^ C2. The explicit expressions for the heat kernel coefficients up 
to 02 we take from Ref. [|19|, where a compact cylinder with unequal velocities of light inside 
and outside was considered. When c\ = C2 the final expressions for these coefficients are 
drastically simplified 



ao = 0, ai/2 = 0, ai = 0, 



03/2 



3f 



(47r)3/2 64 i?' 



02 = 0. 



(5.13) 



The zeta function obtained for given boundary conditions in Ref. |T8| gives 

-2 45 



Q5/2 

(4^)3/2 



e 



i?3 8192' 



0, J 



3,4,5,.... 



(5.14) 



The heat kernel coefficients (|5.13|) and ( |5.14D lead to the following high temperature behavior 
of the internal energy in the problem at hand 



U{T) 



3^2t 
64 i? 



5 C'^^ \ ,^,rr. S^ 



512 



(5.15) 



The corresponding zeta determinant is calculated in Appendix B 

C'(O) = ^ fo.20699 + — In — V 
Now we can write the high temperature asymptotics for free energy 



(5.16) 



F(T) = -e\ 



0.10350 + Am ^ + ^^4^' 

64 2nc 65536 B?T^ 



+ o(r 



-3a 



(5.17) 



and for entropy 



S{T) = ^ 



15 c^fi^ 

0.10350 + — ( 1 + In ^ ) ^^ 

64 V 2hcJ 65536 T2i?2 



3 /. , RT\ 



+ 0{T- 



(5.18) 



Putting in these equations ^^ = 1 we arrive at the double scattering approximation for a 
perfectly conducting cylindrical shell (see Eq. ( |5.10|) ). A slight distinction between the linear 
in T terms in Eq. ( 5.10| ) and Eq. ( 5.17| ) is due to a finite error inherent in the numerical 
methods employed in both the approaches. 

In the case of a pure dielectric cylinder (/ii = /i2 = 1, si 7^ £2) the first four heat kernel 
coefficients are different from zero even in the dilute approximation ||19| (small difference 
between the velocities of light inside and outside the cylinder) 



oo = T- (ci - C2) H 5 — (ci - C2) , ai/2 = 1 — (ci - C2) , 

C2 C2 C2 

14 7r , ,2 3 7r3/2 



-2 '-''^2 -LUJLU2 



ai = ^ (ci - C2) - ^^ (ci - C2) , 03/2 = Y^^^ (ci - C2y 



_ «5/2 _ 857 (Ci-C2)^ 

"^-°' (4^ - 61440 m ■ ^^-^^^ 

It should be noted that the coefficient 02 vanishes only in the (ci — C2) ^-approximation. As 
a matter of fact 02 contains nonvanishing (ci — C2)^-terms and those of higher order [19 



Therefore the zeta regularization provides a finite answer for the vacuum energy of a pure 
dielectric cylinder only in the (ci — 02)^- approximation even at zero temperature. 

The contribution to the asymptotic expansions of the ffist three heat kernel coefficients 
should be involved into the relevant phenomenological parameters in the general expression 
of the classical energy of a dielectric cylinder (in the same way as it has been done for a pure 
dielectric ball). By making use of the coefficients 03/2 and 05/2 we get the high temperature 
asymptotics of the internal energy in the problem at hand 

U{T) = An' — - (1 - -^^^] + 0{T-') . (5.20) 

^ ' 128 i? V 17280 T2i?2y ^ ^ ^ ' 

where An = ni — 722 — (c2 — ci)/c. 

In view of a considerable technical difficulties we shall not calculate the zeta function 
determinant for a pure dielectric cylinder. We recover the respective asymptotics of free 
energy by integrating the thermodynamic relation (|2.11|) and of entropy by using the relation 



( p.l2| ). Pursuing this way we introduce a new constant of integration a that remaines 



undetermined in our consideration 



^(^' = -^"^lis h '■> f + ^^^ 1 + O(T-) , (5.21) 



857 c2;^2^ 


+ C(T-2 


34560 T2i?2^ 


857 c'h'\ 


+ C(T-3 


3A5Q0Tm^j 



S{T) = An^^ (l + « + ln|^-:;^^;^^)+ 0(T-) • (5.22) 



VI. CONCLUSIONS 

In this paper we have demonstrated efficiency and universality of the high temperature 
expansions in terms of the heat kernel coefficients for the Casimir problems with spherical 
and cylindrical symmetries. All the known results in this field are reproduced in a uniform 
approach and in addition a few new asymptotics are derived (for a compact ball with ci = C2 
and for a pure dielectric infinite cylinder). 

As the next step in the development of this approach one can try to retain the terms 
exponentially decreasing when T — > 00. These corrections are well known, for example, for 
thermodynamic functions of electromagnetic field in the presence of perfectly conducting 
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parallel plates [0,0 (see also Eq. ( [4.23|) ). In order to reveal such terms, first of all the 



exponentially decreasing corrections should be retained in the asymptotic expansion 
for the heat kernel. 

It is worth noting that in the framework of the method employed the high temperature 
asymptotics can also be constructed in the problems when the zeta regularization does not 
provide a finite value of the vacuum energy at zero temperature, i.e. when the heat kernel 
coefficient 02 does not vanish. 

In Ref. ||13[ it was argued that in the high temperature limit the behavior of the Casimir 
thermodynamic quantities should be the following. In the case of disjoint boundary pieces 
the free energy tends to minus infinity, the entropy approaches a constant, and the internal 
energy vanishes. Contributions to the Casimir thermodynamic quantities due to each indi- 
vidual connected component of the boundary exhibits logarithmic deviations in temperature 
from the behavior just described. In our consideration we were obviously dealing with an 
individual connected component of the boundary (a sphere or cylinder). Our results corrob- 
orate the relevant conclusions of Ref. [|1^ concerning the free energy and entropy. However 



the internal energy in our calculations tends to infinity like T instead to vanish, this increase 
being caused by the respective logarithmic terms in the high temperature asymptotics of 
free energy. 
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APPENDIX A: ZETA FUNCTION DETERMINANTS FOR 
ELECTROMAGNETIC FIELD SUBJECTED TO SPHERICALLY SYMMETRIC 

BOUNDARY CONDITIONS 

1. A perfectly conducting sphere 

First we calculate C'(0) (zeta determinant) for electromagnetic field in the background 
of a perfectly conducting sphere. We proceed from the following representation for this zeta 
function |jl8| 



C(^)=(-) '-^^^^^Y.i^l + 1) r dyy-'^^Hl-aKy)]. (Al) 



where 



y 



ai{y) = -^[yh{y)KMl v = 1 + 1/2. (A2) 
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The analytic continuation of Eq. ( [AID to the region 3 s < is performed by adding and 
subtracting from the integrand its uniform asymptotics at large v 



t'iz) 



ajiv z) 



4z/' 



2 ' 



t{z) 



VTT^ 



As a result we obtain 



C(s) = Z(s)+Cas(s), 



where 



Us) - [-) 






Zis) = {-j "-^^^^Y.^ 



t i- |ln[l - afiuz)] 



1 = 1 

l-2s / j^ „l-2s .8 



Att 



E 
1=1 



V 



-2^* dz 

dzz'-'^'fiz 



4z/2 ri + z2)3 



(A3) 



(A4) 



(A5) 



2s 



K?)' ^(i+^)(2+^) 





l+2s 



1)Cr(1 + 2s)-2 



l+2s 



(A6) 



When calculating C'(0) o'^s can put in Eq. ( [A 51 ) s = everywhere except for sin(7rs), 
the latter function being substituted simply by vrs. In view of this the integral in Eq. (|A5|) 
is evaluated easy if one takes into account the limits 



limaf(z/z) 

2^0 



Tiu) 



lim 



t\z) 1 



\2T{u+l)) 4z/2' z^o 4z/2 4z/2 
and the asymptotics ( |A3| ) at large z. As a result we obtain 



(A7) 



Z'(0) = -2E^ 



In ( 1 - -—,\ + -^ 



(A8) 



Differentiation of Eq. (|A6|) with respect to s at the point s = gives 

C(o) 

In order to calculate the sum over / in Eq. (|A8|) we consider an auxiliary sum 

a2 



5 1 7 

- + -lni? + ln2 + - 

8 2 2 



(A9) 



1=1 



In 1- 



4z/2 / 4 z/2 



^(0) = 0, 5(1) = Z'(0), 



(AlO) 



where a is a parameter. Derivative of this sum with respect to a can be rewritten in the 
form 



S\a) 



E 
1=1 



+ 



1 + 1/2 l + {l + a)/2 1 + 1/2 l + {l-a)/2 



(All) 



The summation in Eq. ( [A 11] ) can be done by making use of the following relations [^ 
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1 



1 



^^\y + k x + k^ 



1 1 

X y 



+ ^(x) -tp{y), 



^(x+ 1) =^(x) + -, '0 f -j = -7 



- 2 In 2, 



(A12) 



where iIj{x) is the digamma function (the Euler ip function): iIj{x) = {d/dx)\nr{x). This 
gives 



5'(a) = a(2-7-21n2)-- 



, , 3 a\ /3 a 

■^'li + iJ + ^U-i 



(A13) 



Now we integrate the both sides of Eq. ([A13|) over a from to 1 by making use of 'Maple' 

5(1) = Z'{0) = ^-| + ^ln2 + 6 C^(-l) . (AM) 



From Eqs. (^, (^, and (^T|) it follows that 



C'(0) = ^-| + ^ln2 + 6a-i; 



'- + -In— + ln2 + - 
8 2c 2 



-i + ^ln2 + 6Ck(-l) + ^ln- 
8 2 c 

1 p 

0.38429 + -In-. 

2 c 



(A15) 



2. A material ball w^ith ci = C2 

The same technique can be used for calculating the zeta function determinant in the 
case of equal velocities of light inside and outside the material ball (see Section IV.B). The 



complete zeta function in this problem has the form |]T8 



2s 



ny sin(7rs^ °° 
c J vr 



y (2/ + 1) / dyy~'^ — ln[l - ^Vf (y)], (A16) 

t^i Jo dy 



where (Ji{y) is defined in Eq. ( |A2D and the parameter ^^ was introduced in Eq. ( [4.15^ . 
Adding and subtracting under the integral sign in Eq. ([A16|) the uniform asymptotics of the 
integrand at large v we get 



,, , /i?\^"sin(7rs)^ ^_^^ r°° dz d f, ,^ .3 2/ m , '^^ 1 1 



+eus) 



(A17) 



where the function Cas('S) was introduced in Eq. ( |A6| ). Proceeding in the same way as in the 
previous subsection we obtain for the derivative of the function ({s) at the point s = 



c'(o) = 5(a+rCs(o) 



(AI8) 
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where the function 5'(^^) is defined in Eq. (|A10|) . For small values of the argument ^^ we 
deduce from Eq. ( |A10|) 



e 



4 CO 



e 



(A19) 



Therefore restricting ourselves to the first order of ^^ we arrive at the final result 

C'(0) = e^C(0)=e^(-^ + ^lnf + ln2 + |). (A20) 



3. A pure dielectric ball 

A material ball with arbitrary velocities of light inside and outside treated in Section 



IV. C proves to be a more complicated problem. In notations of Ref. ||21[ the relevant zeta 
function takes the form ({s) = (-i{s) + Cil-s), where 



2i?2 



"-^ Tis + l)Ti-s)t, 



d 



1 



^z/i-^M dkk--''^\n^,^i{uk), u = l + ^, p = ±l (A21) 
Jo dk 2 



with 

^p,i{^^) = / p .N [x''s'i{j^h) ei{uk2) - si{uki) e'i{i>k2)\ , fci = fc/ci, k2 = k/c2, (A22) 



si{y) = \hrlu{y), ei{y) = J — K{y). 

V Z V TT 



The parameter x = v (^i/^2)/(^2/^i) corresponds to ^ in Ref. pl[ . 

The analytic continuation of the zeta function at hand to the region 9fJ s > is performed 
by adding and subtracting in Eq. ( [A22| ) several terms of its asymptotic expansion 



Ap,K^^)~ E 



D 



n,p 



(A23) 



n=-l,0,l 



where 



D-i = r]{ki) - 7]{k2) - (ki - fca), vi^) = VT+z^ + In 



1 + VTT 



Dn = ln 



X^ Ci t2 + C2 ti 



\/ClC2tlt2(X^+l) 



1, 2. 



1 + kf 



(A24) 



For our purpose it is sufficient to consider four leading terms of the asymptotic expan- 
sion ( [A23| ), n = —1, 0, 1, 2. The functions Di and D2 are given in Ref. |21|. Proceeding in 
this way we represent the zeta function ( |A21| ) as follows 
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(A25) 



where Ch is the Hurwitz zeta function. Taking the derivative of the zeta function ( |A25|) at 
the point s = with allowance for the behavior of Di{k) at k = and A; = oo we obtain 



°° / 1 

c;(o) = -2e(^ + ^ 



^l 2i/x''ci + C2y 2z/x''ci + C2 8z/2 lx''ci + C2 



+2 li„£i + iii„ 

1 4 ci 24 



X^ Ci + C2 



+ 



1 X^ Ci - C2 

2 x^ ci + C2 



1 



;2-lni?-7-2 ln2) 



X^ Ci - C2 



CO d 

dk \n.k-— D2 

dk 



(A26) 



In order to calculate the sum over / in Eq. (|A26 ) we consider an auxiliary sum 

5i(0) = 0, 



1=1 



\ f, b\ b b^ 

In 1 + :—-;— + 



2uJ 2z/ ■ 8i/2 

with b being a parameter. The derivative of this sum can be cast to the form 

1 1 



(A27) 



b °° 

^ 1=1 1 



1 + 1/2 / + (& + l)/2 
Taking into account Eqs. ( |A12|) we obtain 

5'(&) = ^(-2 + 7 + 21n2) + ^V^[^ + ^ 



(A28) 



(A29) 



The integration of Eq. ( [A29| ) over b from to 6 = {x^ Ci — C2)/{x'' Ci + C2) gives the sum 
entering Eq. (|A26|) 



5i(&) = j(-2+7 + 2 1n2) + 6CH(o,^ + ^ 



-2 



'^"(-'■l + ^) + '^"(-4 + ^^ 



+ 2 



ch(-i.|)h-c;,(-i,5 



(A30) 



Substitution of Eq. ( [A30|) into Eq. (^) gives 

C(0) = -^(-2 + 7 + 2 ln2)-K^(o,^ + ^ 
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b b^ 
+ ---(2-lni?-7-21n2) 



C2 



X'' Cl + C2 



(A31) 



In the case of nonmagnetic media (yUi = fi2 = 1) the right-hand side of Eq. ( |A31| ) is 
shghtly simphfied. Assuming that we are deahng with a dilute dielectric ball we can expand 
C'(0) in powers of the difference (ci — C2), where ci = l/i/ei, C2 = 1/a/^- As a result we 
get 



C'(0) = C;=-i(0) + C;=i(0) = ^(-^ + ln| + ln4 + 7)^^^ 



+ Oi{c,-C2f). (A32) 



APPENDIX B: ZETA FUNCTION DETERMINANTS FOR 
ELECTROMAGNETIC FIELD WITH CYLINDRICALLY SYMMETRIC 

BOUNDARY CONDITIONS 

1. A perfectly conducting cylindrical shell 



A complete spectral zeta function in the problem at hand is defined by the following 
expression 



C(^) 



20Fc-r(.)r(3/2-.)io dyy'-^'^yHi-f^liy)] 



K 



2s-l 



0Fc2^r(s)r(3/2-s) 



Y: n'-'^ / dyy'-'^— ln[l - f^l{ny)], (Bl) 

n=i ^0 dy 



where 



f^niy) = y^[In{y)Kn{y)]. 

The first term on the right hand side of Eq. ( plD is an analytic function of the complex 
variable s in the strip —1/2 < ^s < 1/2. Therefore there is no need in analytic continuation 
of this expression when calculating ('{0). As regard to the second term in Eq. ( pl|) its 
analytic continuation to the region ^s < can be accomplished in a standard way. We add 
and subtract here the uniform asymptotics of the integrand when n tends to infinity 



ln[l - fJ'liny)] 
As a result we obtain 
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^ 1-2. r ^^ ^ l^ ri 2. M^l/'^'l 



0rc2^r(s) r(3/2 - s 

i?2«-i r(i/2 + s) 



Keeping in mind the behavior of the gamma function at the origin T{s) ~ s~^ one can easily 
find the derivative of ({s) at the point s = 



2 
4 



Unhke the spherically symmetric boundaries, the integration is not removed in the formula 
obtained for C'(0). Therefore the first two terms in Eq. (B4) can be calculated only numer- 
ically 

^ r ^ ^ I, 2! M 0.53490 ,^^, 

-— / dyHl-iil{y)] = ^. (B5) 

Tr/t JO H 

Applying the FORTRAN subroutine that approximates the Bessel functions by Chebyshev's 
polynomials we evaluate the first 30 terms in the sum in Eq. 



2 ^ /-oo r, ^^ 2^ ^^ yH^\ 0.00554 



Finally gathering together all these results we have 



045847 ^,_^fi 
^^ ' R ?,2R 2c ^ ' 



2. A compact infinite cylinder with ci = C2 

Now we turn to a compact cylinder placed into unbounded medium such that the velocity 
of light is uniform on the lateral surface of the cylinder. Proceeding as in the case of a 
cylindrical shell we start with the expression for the relevant spectral zeta function 

p2s-l oo oo J 



with parameter ^ determined in Eq. ( ^.15| ). In the linear approximation with respect to ^ 



2 



Eq. ( pSf) assumes the form 
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The analytic continuation to the region 3?s < is needed only for the second term in Eq. 
M). Adding and subtracting here the uniform asymptotics of the integrand for large n 
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we obtain 
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The first two terms in Eq. ( Bll ) can again be calculated only numerically 
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The final result reads 


C'(o) = 1 


0.28428 - 0.00640 + — 


e 


( n on/^n 
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